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1. Introduction. The model and some relations between its parameters 

String breaking at zero and finite temperatures is known as one of the important open 
problems in QCD (see e.g. refs. [Q, ^, ^ |5|, 0] for recent developments). The aim 
of this paper is to proceed further with the calculation of the string-breaking distance, 
along with the lines of ref. |^] (Sections 2-6), as well as using some alternative model 
(Section 7). In ref. dual Abrikosov-Nielsen-Olesen strings have been used to model flux 
tubes of the chromoelectric field in real QCD. Such strings Q are solutions to the classical 
equations of motion in the 4d dual Abelian Higgs model, whose Euclidean Lagrangian reads 
£ = IF^, + \D^^\'' + |(|99|2 - z;2)2. Here F^, = d^B, - d,B^, D^ip = {d^ - ign,B^)ip, B^ 
is the dual gauge field, if is the complex-valued dual-Higgs field, and gm is the magnetic 
coupling constant, related to the electric one, g, as gm = "^T^/g- The masses of the dual 
vector boson and the dual Higgs boson are my = V^gmV and rriH = \/2Af , respectively. 

We will consider this model in the London limit ^, L = In ^ 1, where all the 
results can be obtained analytically. The Bogomol'nyi case, L = 0, will be considered 
elsewhere. The electric field of a straight-line dual Abrikosov-Nielsen-Olesen string reads 

^The ratio is called Landau-Ginzburg parameter. 
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E{r) = Y^Ko{'mvr), where r = |x_l| with x_l = {xi,X2) denoting the direction transverse 
to the string, and from now on Ky^s stand for MacDonald functions. The field averaged 
over the string cross section, (E) = ^ f (PrE{r) obeys the relation g {E) = Aa/L. Here 
S = TTrriy^ is the area of the cross section of the string, and the string tension has the form 
a = 2'Kv'^L. In order to establish a correspondence to QCD, we will need to express all the 
results in terms of a and g. Two formulae useful for this are niy = and S = 



2. The string-breaking distance with the neglection of dispersion of E{r) 

The rate of the pair production reads w = ^Imr[Aj], where the one-loop effective action 
of a scalar quark has the form: 



dT 



e-'-'^x 



/ 



Px_Ll?X|j exp 



dT -X 



+ - ^■E{y.^{T))eijXiXj 



(2.1) 



Here, xy = {x^,x/i), i,j = 3,4, and the field of the string is Ai = —^eijXjE['K±). Note 
that, in this paper, we will study the case of scalar quarks only (except of some com- 
ment at the end of the next Section), since the longitudinal part of the spin factor. 



/ dr (^tpiipi - igipiil^jFijir) 



where Fs4^{t) = — -F43(t) = E{x±{t)), cannot 



/ Vil^i exp 

be calculated exactly as long as Fij is r-dependent. To evaluate the path integral ( |2.lD , we 
will use the requirement that the mass of the produced pair, m, should be much larger than 
my, i.e. m ^ i'^ order that the field of a string can be considered as a constant 

one. The "largeness" of m means that characteristic proper times are "small", T < 
that enables us to evaluate the path integral semiclassically and to compute further the 
leading quantum correction using the Feynman variational method Q. Furthermore, we 
naturally assume that not only the Compton wavelength of a produced pair, m~^, is much 
smaller than the range of the field localization, m^^, but also that the characteristic pair 
trajectories are small compared to voy^ . Since, in the Euclidean space-time, pair trajecto- 
ries are circles of the radius Kp = ^, the condition of smallness of the pair trajectory, 
Rp <C m-y^, yields m <C 2g^ Both conditions. 



2 /vro" / a 

-W— <m<2gW— , 
g \ L V 



(2.2) 



are compatible to each other at g ^ 1. 



^This can be seen either by solving the respective Euler-Lagrange equation |lO|, or simply by noticing 

where <1> is the flux of the electric fleld through 

the contour of a pair, $ = 3 (£) ■ ttR;. Therefore, a pair is identified with a quark which moves along a 
circle of the radius R„. 



that exp ( — ] in the Schwinger formula should be e 
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Then, due to the smallness of a pair trajectory, E{x±{t)) can be replaced by its value 
averaged along the trajectory: 



1 1 

J dTE{xj_{T))XiXj -J^ij ■ ^ j dTE{x±{T)), 



(2.3) 



where Sjj = J drxiXj is the (i, j)-th component of the so-called tensor area of the trajec- 



tory. In the leading small-T approximation, we obtain the classical result for j 'Dx± = 

J d'^x±{0) J Px_l(t) in eq. (|2l|) : 

xx(0)=x_l(T) 



d'^x± exp 



(2.4) 



Thus, 



iVciV/ f dT^_^.^ f ^^^^ I _1 I ^^^j I I ^2^^^^p 



47r J r2 





E{x±)eijTl^J 







(2.5) 



Neglecting for this Section the dispersion of the field -E(x_|_), we have 



d x± exp 



S ( exp 



E(x±)eijT,ij 



5* exp 



{E) eijEi 



(2.6) 

In this approximation, we therefore arrive at the Euler-Heisenberg Lagrangian in the con- 
stant field Ai = —^CijXj (E), 



N,Nf fdT^_^2T g{E) 



sin(<7 {E) T) ■ 



(2.7) 



and recover for w the 4d Schwinger result in the bosonic case: 



w = NcNf 



- (-1 



■ exp 



TTkm^ 



(2^)3 fc2 g^E) 



(2.8) 



We can further express the inequality T < in terms of the parameters of our model. 
Namely, since at least the first pole in the imaginary part of the Euler-Heisenberg La- 
grangian should give its contribution, T may not be arbitrarily small, but the following in- 
equality should rather hold: T > = The condition ^ > T then yields m < 2^ 



that also coincides with the condition for w not to be exponentially small. This new con- 
straint is stronger than the above-obtained one, expressed by the right inequality of (^), 
since g ^ 1 is now absent. The new constraint can be viewed as an upper boundary on L: 



'Up to the factor NcNf absent in the electromagnetic case. 
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4 a a , ^ 

L < ^ ~ 1.27^. (2.9) 

Setting a = (440 MeV)^ and a typical hadronic mass m = 200 MeV, we get an estimate 
L < 6.2, that still leaves a window for L ^ 1. 

Approximating the whole sum ( p.8D by its first term (equal to the density of produced 
pairs), we have 

2NcNf /cr\2 / vrm^LX 
The respective string-breaking distance r = has the form 

g^NcNfa \ 8 a J 
where, due to inequality (^), exp (^^) < \/e — 1-65. 



3. Finite-temperature generalizations 

Setting for the mass of a produced pair the vr-meson mass, we can further extend the 
analysis of the previous Section to the case where the temperature is close to the critical 
one. This can be done by using the formulae |^, |lT[ m? = m^t^'"^^, a = aot^'^"^, where 
t = 1 — ^ is the reduced temperature ^, and do ~ (440 MeV)^. We obtain: 



w 



^(l)^^O (3.1) 



that establishes the law by which r grows at t ^ 0. In the same limit t ^ 0, the condition 

9 fT 

2 V TO-' 



mr_|_ » 1 with the temperature-dependent r±, r± = f \/ — , leads to the following boundary 



1 / 1 / — 

on t from below: -r -^1 / — ^ t. On the other hand, the condition of smallness of the 

L yrun y as J ' 

pair trajectory, expressed by the right inequality of (2^), sets a boundary on t from above: 
/ 2 \ 1-82 

t TTTUT ( ■ These two conditions imposed on t are apparently compatible to 



each other at sufficiently large L and/or g. Note also that the condition ( |2.9| ) becomes 
softer as one approaches the critical point, since = -^t^^-^^ 00. 

Let us now evaluate the string-breaking distance at relatively low temperatures, namely 
at T = ©(/tt), where the following formula holds p 



mi{T)^ml{l + ^]. (3.2) 



*The fact that the temperature and the proper time are denoted by the same letter "T" should not lead 
to reader's confusions. 
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We can further approximate my(T) by my(0), since rriy is the vacuum correlation length, 
whose QCD-analogue at such temperatures can be considered as temperature-independent 
11 

Next, we need to know the temperature dependence of g {E). Since we are exploring 
the region of temperatures smaller than the temperature of dimensional reduction, the sum 
over Matsubara frequencies appears: 

g{E) (T) = / <fz 2 Ao iyz^TWM^) , 

n=~oo 

where /? = 1/T. To carry out the integral, it is convenient to transform the sum as follows: 



+00 



J2 Uz^ + {mvPny 



1 



+00 



2 E 

n=— 00 



dt 

-exp 



---t{z^ + {mvPnf) 



2mv 



E 



dt 



exp 



_^2^ _ ]_+ (Wn/mv)' 



4t 



my 



(3.3) 



where lOn = 27rTn. The integration over d^z then immediately yields 



g {E) (T) = ^ 



[1 + (a;„/mv)2] 



213/2 • 



(3.4) 



To study the zero-temperature limit of this expression, one should perform the inverse 
transformation of the sum: 



+00 



[1 + {urjmvff^ 



^ y +00 



{my(3n 



i2n 



{mvP) ^ |n|i^i(mv/3|n|). 



vr 



At small T's of interest, the sum here is apparently dominated by the zeroth mode. More- 
over, expanding the function Ki in power series, one can check that no finite corrections 
exist to the r.h.s. of the formula |n|iri(my/3|n|) ^-^-^ T/my. Therefore, the leading small- 
T (physically, at T ^ f^^) correction to g {E) (T) stems from the terms in the sum with 
|n| = 1, 2Ki{mYf3). This correction is therefore exponentially small, as confirmed by the 
following final expression: 



^In QCD, the (magnetic) vacuum correlation length becomes definitely temperature-dependent at tem- 
peratures larger than the temperature of dimensional reduction (that is of the order of 2Tc), where it is 
1 

9^(T)T- 



proportional to ^ 
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(3.5) 



where o"o = 2'kv'^L. As for the dependence cr{T), it is derived in Appendix B. 



Equations ( |3.2| ) and ( |3.5| ), being substituted into the formula w ~ N^Nf (^^yi e 
determine the temperature dependence of tu. Apparently, the correction produced by 
eq. ( p. 51 ) is exponentially small with respect to the temperature dependence appearing by 
means of eq. ( |3.2| ). Therefore, w decreases at T increasing from zero to the temperatures 
of the order of 7,^. The corresponding increase of f parallels the same phenomenon we have 
found at T ^ Tc- 

It is finally worth making the following comment. In the case where the dispersion 
of E(r) is neglected (that we are discussing here), the spin factor can be evaluated, and 
it is natural to address the issue of how strongly the antiperiodic boundary conditions for 
spin-^ quarks can affect the obtained result. For such quarks, one should substitute in 
eq- (£D sin{a}E)T) ^ - tan(g(i;)r) ' that yields, instead of eq. (P), 



w 



47r3 



k=l 



exp 



nkm 
'J{E) 



(3.6) 



It has been shown in ref. |14] that antiperiodic boundary conditions for fermions can 
be taken into account upon the multiplication of the zero-temperature one-loop effec- 



tive action by the factor 



, where (only for this formula in this 



1 + 2 X;(-l)"e- — 

n=l 

Section) T stands for the proper time. In course of taking Imr[Aj], this factor trans- 
forms to another one, by which a k-th term of the series in eq. (|3.6| ) should be multiplied: 

9(E) 



1 + 2 ^(-l)"exp 

n=l 



9{E)l3^n^ 



For relevant /c's, which are k < we obtain the 



following relevant n's: n < 2T 



irk 
9(E) 



2T 



that, at T < /tt, are smaller than 2. Therefore, 



the factor produced by the antiperiodic boundary conditions for spin-i quarks reduces to 



1 - 2 exp 



9{E)f3^ 
ink 



Since, for the above-mentioned relevant /c's. 



9{E)P^ 
ink 



> 



~2~ I ' 



T < ^ the obtained correction falls off as the tail of the Gaussian distribution ^. 

v2 



4. Significance of corrections due to the dispersion of E{r) 



Let us now take into account the second cumulant in eq. ( |2.6D [i.e. the dispersion of E{r)] 
and also address the issue of convergence of the cumulant expansion. We have 



exp 



exp 



{EY 



"The inequality T < stems from the fact that the width of the Gaussian distribution e 



V2 
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vrc 



I dfe^p , (4.1) 



where 

^ {E^) - {Ef^ . (4.2) 

The last formula in ( |4.lD means that, instead of the constant field (E), we are now dealing 
with a shifted (but still space-independent) one E = {E) — where the field / should be 
eventually averaged over. The respective pair-production rate reads 



w 



+00 _ , ^ 

— ^ / dfe'^ ■ -f— V > ^^-r4 — exp ^ . 4.3) 

-00 '^—^ 



The parameter of the cumulant expansion, i.e. the ratio of the absolute value of the second 
cumulant to that of the first one, is 



Here, 

(E^) - {Ef , , 

d = ^ — ' ; ' 4.5 

{Ef 

is a measure of dispersion of the field E{r). It is calculated in Appendix A and reads 
d = ^ — 1. The condition of convergence of the cumulant expansion, i.e. the demand that 
the parameter (^^) is smaller than unity, then yields yet another upper boundary on L: 

4 a a , , 

L<— — -2.23— , (4.6) 



that is, however, weaker than the condition ( |2.9D . As well as eq. ( p. 9]) , this new condition 
on L becomes softer at T — > Tc, since = -^t^^'^^ — > 00. In another words, at a 
fixed L, obeying condition ([4.6|), cumulant expansion converges better at T — > Tc, since its 
parameter vanishes in this limit as 0{t^'^^). 

To perform the average over / in eq. ( f4.3D , note that the dispersion (i.e. the width) 
of the Gaussian distribution, ^ is Therefore, characteristic /'s obey the estimate 

\f\ < V2c = f^(£;2) - {Ef, and < Vd ^ Vol7 ~ 0.76. Therefore, although 

with the accuracy of only 76% (cf. Section 5, where this problem will be avoided), we can 
approximately write 



When this expansion is substituted into eq. ( [4.31) , the /-integration can already be per- 
formed and yields [cf. eq. (^^)] 
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w 



{9{E)f 
(27r)3 



E 

fc=i 



2Tim?k 



exp 



■Km k 



irm k 



(4.; 



Approximating again the whole sum by the first term only, we see that the obtained 

2 

correction is small, provided <^ 1. This is precisely the condition of convergence of 

the cumulant expansion, eq. (4^), with "<" replaced by "<C". In particular, the correction 
produced by the second cumulant becomes vanishingly smaller than the leading term at 
T — > Tc. We also see that the obtained correction increases w, thus diminishing f. Due to 
condition (|2.9|) , f diminishes only by a factor of the order of unity. 



5. Accounting for quantum effects by the Feynman variational method 



In this Section, we will evaluate the leading quantum correction to eq. (|4.8| ). It can be 
obtained upon the small-T analysis of the path integral over x±(t) in eq. ( |2.lD with the 
approximation ( |2.3| ) adopted. Such an integral can naturally be evaluated using the Feyn- 
man variational method In 2d-case and in our notations, it looks as follows. We need 
to evaluate the path integral 



Z{T) = j (fx{id) j Vx{T)e-^, 



x(0)=x(T) 

where S = J dr + U{x.)j , U{x) = ^E{x)eijY,ij. The classical expression for this 
^ ^ 

integral, J (Pxe~^^^^\ given by eq. (2^), corresponds to T's, which are so small that 
the trajectory does not deviate from its initial point x(0). Let us further introduce the 

T 

coordinate describing the position of the trajectory xq = y / (iTx(T), the trial action 



Sq = f dr^ + TW{xq), and the respective trial partition function 





d^XQ 



Pxe 



-So 



fixed Xi 







d^xo 



Vye 



1 



AirT 



(5.1) 



y(0)=y(T)=0 

where y(T) = x(t) — x(0). Here TV is a trial function, which should be determined upon 
the minimization of the expression Fq -\- ^ {S — ^0)^^^ (which approximates the true free 
energy F of the system from the above, F < Fq + ^ {S — 'So)^^), where F = In^(T), 
Fq = — yln2^o(r)- It is further possible to demonstrate that the averaged difference of 
actions can be written as 



1 



0/-S0 



So 
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For the first of tlie two averages here one has ({7(x(0)))^^ = / d?ylC{y)e ^^(y), where 

/C(y) ^ I V^Ui^m exp I - y dr^ j 5{y - xq) = 
x(o)=xm \ / 



(2^ 



-iky 



Px exp 



L 



x(0)=x(T) 

and U{q) = f d'^xe^^^^U{x) is the Fourier image of C^(x). Using the formula 

T 

Px exp 

x(0)=x{T) 



X 



dr { — + ifx 



vr 



-(5 I / drf I exp 



T T 2 / \ 

i 1 dr 1 dr'lr - r'|f(r)f(r') + ^ E / ^^^/^ 
\0 / 



with f = y + q(^(r), one obtains /C(y) = -p^^^ / c?^xC/(x)e y)^. Thus, according to 
eq. Q, 

^ _ /ciV-^^(y) [4^r/C(y) - Ty(y)] 

Note that the potential smeared over the Gaussian distribution, 47rr/C(y), goes over to the 
unsmeared potential U (y) in the limit T ^ 0, as it should be. 

The variational equation 6 [Fq + y (5 — ^o)^^] = then yields W{y) = 47rT/C(y) as 
the best choice of W{y) . Accordingly, the new expression, which accounts for quantum 
effects, 



AirT 



d^ye 



4ttT 



d yexp 



j d^xU{x)e-^ 



(x-y)2 



(5.2) 



is a better approximation to Z(T) than its purely classical counterpart J d'^xe~^'^^^^ 
(recoverable in the limit T ^ 0), which was used before. 

The obtained result (^^) prescribes to replace E{-k±) in eq. (|2.6|) by 



ttT 



d^yE{y)e- 



In particular, {£) = {E), i.e., at the level of the first cumulant, the variational method 
yields the same result as the clasical approximation. We further obtain 



^With this choice of W(y), {S - So)^^ = 0. 
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This expression replaces {E'^) in eq. (|4.1| ); in particular, (<?^) — — ^ (-^^)- Thus, within the 
variational approach, the constant (^^) becomes replaced by Cvar = ( ('^^) ~ 



8 

apparent difference of Cvar from c is that the former is T-dependent, whereas the latter is 
not ^. The one-loop effective action reads 

oo +00 _ 



32^5/2 J 2-2 J J sm{gET) 







where "~" means "in the bilocal approximation to the cumulant expansion", and again 



2/ 
a ' 

Let us further calculate the integral entering Cy 



E^{E)-f 



[ (fxd^ye-^^''-y^'Ko{mv\^\)Ko{mv\y\) = ^ [ d^zd^z'e ^ A Koiz)Koiz'), 

(5.3) 

where z = |z|, z' = |z'|. This can be done approximately, by using the fact that T < ^ 
i.e. TfT-^ ^ 1, that allows us to Taylor expand Kq{z') around z. This calculation, 
whose details are given in Appendix A, leads to the following result: 



3 

where the variable y = gET acquires the values vr/c when one is taking Imr[^j]. We can 
further use the approximation ( ^ ) with the term neglected, since it would otherwise 
lead to the excess of accuracy, making the /-integral non-Gaussian. Denoting this approx- 
imation by "~" and substituting the above-mentioned values of y, we obtain the following 
modification of eq. (4.2): 



9' 



This finally results in the following change of the dispersion parameter (^]^), entering 
eq. (|4.8[), that makes this parameter ^-dependent: 



d ^ dk 



{Ef 



This formula describes the quantum correction to eq. (4^). We, however, see that this 
correction is small. Indeed, due to the main exponential factor in eq. (|4.8| ), we have for 
relevant /c's: 



*In particular, this means that the /-integral in eq. ( |4.l[ ) should now be considered only inside the 
T- integral (and not vice versa). 
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irk TT g (E) a 1 

where the last inequaHty stems from the condition m ^ mv [expressed by the left inequality 

of (PD]. 

6. Accounting for the space dependence of E{r) without the cumulant 
expansion 

In this Section, we will evaluate w in an alternative way, namely completely without the 
use of cumulant expansion in eq. (l2.6|) . The necessity of this calculation is, firstly, because 
cumulant expansion has been shown to converge only provided L is bounded from above 
according to the condition (|4.6|). The new method enables one to replace this constraint 
by some other, weaker, one, and to relax the constraint ( |2.9D . Secondly, the approximation 
used to arrive at eq. ( [4.8D , which was holding with only 76% accuracy, will now be not 
necessary anymore. 

The method of this Section is based on averaging every term in the expansion of 
the exponent in eq. ( |2.6| ) over d^xj,. In the London limit, such an average can be done 
analytically by making use of an explicit form of E{r). Namely, we have 



d x±_ exp 



The last integral approximately equals -^-2^ "n!, that yields 



47r 



E 

y n=0 



4 Qm 



4tt 

1 + ^^.S^ ^ 



dte 



t 1+ 



-22 \ 



The series above converges at • 2nRp < 1, that results in the following boundary on 



L from above [cf. eq. (ffjq)]: 



We further obtain for eq. (2.5): 



16 a a 
L < ^ ~ 5.09^. 



(6.1) 



IT 



^e-™ ^ / Px,i / dte 



s- 



NcNf fdT^_^2j. 



47r2 J r2 




dte 



-t 9E 







sm{gET) ' 



(6.2) 



where E = is a new space-independent electric field. We see that the constraint (^ 
is now removed. Indeed, the condition for the first pole in Imr[Aj] to contribute, vr < gET, 
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together with the inequahty T < yield vr < ^ = . This condition does not produce 

2 

anymore a constraint on L, but merely means that fs obeying the inequality t > ttL^ 
give a dominant contribution to Imr[Aj]. 

The pair-creation rate stemming from eq. ( |6.2[) takes the form 



^ oo 

fc=l 







Clearly, only terms with k < ^^^^-i are relevant in the sum. At T = 0, substituting 
a = (440 MeV)^ and a typical value m = 200 MeV, we obtain k < < 1. Therefore, only 
the first term is relevant, that yields 

^ - ^2^^^^ . (6.3) 

According to the above-obtained constraint ( |6.lD , the argument of the MacDonald function 
in this formula is smaller than 8. For the values of L, at which this argument is still larger 
than unity, i.e. L > 

^5/2^3/4 1^ 

When L obeys simultaneously inequality (|2.9|), the obtained expression can be compared 
to eq. ( 2.10| ). The parametric dependences of these two expression on m, a, and L are 
apparently different from each other. The regime under discussion, 2mA > 1, however, 



does not imply the exponential smallness of (|6.4D , since, due to (1^, 2myjT^ < 4. The 
string-breaking distance, corresponding to eq. (|6.4|), reads 



Its ratio to distance (2.11) is given by the function 1/4 .5/4 * , where the variable 



X = m\l^ ranges between ^ and 2. This is a monotonically decreasing function, which 



therefore acquires its maximum at x = ^ [where approximation (6.4) to eq. ( |6.3| ) starts 



breaking down]. The value of the maximum, ~ 9.10, is, thus, an upper limit for the ratio 
of the two string-breaking distances. The minimal value of this ratio, corresponding to 
X = 2, is only ~ 2.84. 

Instead, at T ^ T^, ^ = ^t^'^^ ^ 0, and we have w (f )^ that differs from 

eq. ( |3.1| ) only by a factor 2. The respective string-breaking distance is, therefore, larger 
only by a factor \/2. 
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7. Considering pairs as holes in the confining peUicle 



In this Section, we will consider an alternative approach to the pair production, based on a 
combination of the formulae on the metastable vacuum decay [l^ with the stochastic 
vacuum model (for a recent review see |]l8| ) . The idea is to consider the produced pairs 
as holes in a 2d confining pellicle, which spans over the contour of an external qq-pair. Such 
a hole is a region where the pellicle is eaten up. Therefore, a hole of a radius R diminishes 
the action of the pellicle by a ■ vri?^, but increases this action by (m + /i) • lirR. Here 
/X is some parameter of dimension [mass], which is a nonperturbative part of a constant 
entering the perimeter law of the small-sized Wilson loop of a produced pair. The critical 
radius (i.e. such a radius, that all holes with R < Rc collapse, while those with R > Rc 
expand and destroy the pellicle) stems therefore from extremization of the action S{R) = 
(m+//)-27r-R— cT-vri?^. This critical radius and the action of a critical hole are -Rc = {m+fi)/a 
and S{Rc) = 7r(m + fi)'^/a. Accordingly, in this approach, the rate of the pair production 
is proportional to e~'^(™+'^)^/'^. The proportionality coefficient [16|, should furthermore 
be multiplied by the factor Nf ^: 



w 



a 

Nf — exp 



(7.1) 

Next, to evaluate both n and a within the same model, it is natural to use the stochas- 
tic vacuum model. It yields the following expression for the Wilson loop (W^(C'))ym = 



c 



YM 



(W^(C^))YM-^trexp 



1^ 
2! 4 



da 



E{C) 



(x) J daxp{x 
S(C) 



YM 



(7.2) 

where S(C) is the surface encircled by the flat contour C. Next, = A'^T"', where 
T°'s stand for the generators of the group SU(A'^r;) in the fundamental representation, 
[r",r^] = if^b^T"", tr TT* = 15"^; the average (. . .)ym is implied with respect to the 
Euchdean Yang-Mills action, \ f d^x{F^^f, where F^^ = d^Al - duA^ + gf^'^A^^Al, a = 



1, 



1; $2,3,/ = exp ig J dz^Ap{z) is a phase factor along the straight line. 



which goes through x' and x. The symbol "~" in eq. (7.2) is implied in the sense of the 
bilocal approximation to the cumulant expansion. This approximation, supported by the 
lattice data |15, 19, |2^], states that, in the Yang-Mills theory, the two-point irreducible 
gauge-invariant correlation function (cumulant) of -F^i^'s dominates over all cumulants of 
higher orders, which are therefore neglected. Finally, the factor 1/2! in eq. ( [7.2p is simply 
due to the cumulant expansion, whereas the factor 1/4 is due to the non-Abelian Stokes' 
theorem. 

^The factor Nc is now absent, since the Wilson loop, describing a produced pair, is a colorless object. 
Instead, the number of Wilson loops, which can potentially be created, is proportional to the number of 
different quark species, i.e. to Nf. 
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The stochastic vacuum model suggests further the following parametrization of the 
two-point cumulant: 



+ 9 [d^l {{x - x')xS,yp -{x- x')p6,yx) + df, {{x - x')p5px - {x - x')x6pp)] x 



((x - x'Y) + 



NcC2 



Mn'^lx - x'\* 



Here C2 



(7.3) 



is the quadratic Casimir operator of the fundamental representation, and 

((-^m^)^)ym 10 

241D(0)+Di(0)| 



J\f is the normalization constant, which in 4d reads M = 24|d(o')+d^(o) ] Inserting this 
parametrization into eq. ( [7.2| ) we obtain (applying Abelian Stokes' theorem to the part 
containing derivatives) : 



9^ 

8N, 



dxn d> dx, 



2 j dapuix) J dap^{x)D {{x - X 
■ S(C) E(C) 

NcC2 1 



+ 



G{{x-x'f) + 



C 



C 



(7.4) 



where G(x^) = J dtDi{t). The one-gluon-exchange contribution to the Wilson loop, 



exp 



Sl^^dxp^dx'^ 



1 



exp 



Sna 



is known [21| to yield the renormaliza- 



c c 

tion of mass of a produced pair. Here a stands for an inverse UV cutoff; a <^ Tg, and 
-'s — 



1 GeV ^ is the correlation length of the vacuum 0, [|, |19|, |2§. The asymptotics 
of the Wilson loop at v^|S(C7)| > Tg, where |S(C7)| is the area of S(C), is {W{C))y^ ~ 
g-CT|s(c)| _ 'pj-^jg asymptotics is obeyed by the Wilson loop of an external qq-paii. Instead, at 
V^|S(C7)| = 0{Tg), i.e. for the Wilson loop of a produced pair, {W{C)) 



YM 



^{m+p)L{C) 



where L{C) is the length of C As follows from eq. ([7.4|), the parametrization ( [7.3| ) 
is chosen in such a way that the function D yields the area law, while the term with 



the function Di contributes to the perimeter law. Namely, the string tension reads |24| 
£7 = J d^xD{ii?), while for the perimeter constant ^ we obtain in a way similar to the 

^ f dCG{C^). A derivation of 



Coulomb [|21| interaction between points lying on C: /Li 
this formula is presented in Appendix A. 



^°Note that the one-gluon-exchange contribution, represented by the 
considered as a part of the function Di |18 



-term, can alternatively be 



Here, we rather consider this perturbative contribution 
separately, so that -Di(O) is a finite quantity. 

^^Below, we will see that the typical size of the Wilson loop of a produced pair is indeed 0{Tg). 
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In what follows, we will adopt the exponential parametrization of the functions D and 
L>i PI, III, 13 D{x'^) = Z)(0)e-l^l/^s, Di{x^) = L'i(0)e-l^l/^f . It yields the following values 



of the string tension and the perimeter constant: a = ^^"^ T^ -D(O), // = |j^T^ 1)1(0). 
Introducing the parameter 7 = L'i(0)/D(0), whose lattice value in full QCD is 0.13 it 0.08 
(see e.g. p^), we can rewrite the obtained results as 



For the critical radius of a hole we have Rn = '""'"^ = H -)_ ^Ta. Substituting again 
a = (440 MeV)^, m = 200 MeV, and = 1 GeV, we get for the worst case of the above- 
quoted lattice value of 7, 7 = 0.21, ^ ~ 1.07 This ratio is certainly 0(1), that justifies 
the use of the perimeter asymptotics for the Wilson loop of a produced pair. 

Comparing the pair production rate, eq. (|7.l|) , with eq. ( 2.10|) , we see that the pa- 



rameter (1 + ^)^ replaces the parameter L/4 of that equation. At o" = (440 MeV)^, 
m = 200 MeV, we have (l + ~ 1.08, whereas, according to eq. (|]9|), ^ < 1.55, that 
is of the same order of magnitude. Finally, the new string-breaking distance, at which the 
potential V{r) = are~'^^ acquires its maximum, is 



r-=^ = ^^exp ^ 1 + ii - . (7.6) 

The argument of the exponent here approximately equals 0.70, that is quite similar to 0.5 
we had as an upper boundary in case of eq. ( |2.11 ). However, contrary to that equation, 



eq. ( [7.61) does not contain the factor yL in the preexponent. 

Notice also that, due to the fact that a = 0{Nj^), the string-breaking distance (|2.1l| ) 
is 0{N^) too Q. For the same reason, as can be seen from the first of eqs. (17.51) , Tg is 
0{Nj^) as well. Therefore, according to the second equation of (^), fi is also 0{N^). This 
means that, as well as eq. ( p.ll| ), the new string-breaking distance ( ^ ) is 0{N^). 



8. Summary 

In this paper, we have considered two approaches to the problem of string breaking in 
QCD: one of these is based on the dual superconductor model of confinement, and the 
other one - on the stochastic vacuum model. In the first approach, which was proposed 
already in ref. [P, the pair-production mechanism is due to the field of the chromoelectric 
flux tube (dual Abrikosov-Nielsen-Olesen string in the London limit of the dual Abelian 
Higgs model). In the second approach, pairs are considered as holes in a pellicle, which 
confines a test gg-pair. 

Within the first approach, we have demonstrated that the result of ref. @, based 
on the Schwinger formula, accounts only for the first term of the cumulant expansion 
in the average ( [4.1| ). In Section 3, we have found the temperature dependences of the 

Notice that Tg itself is smaller (by a factor of the order of 5) than a typical size of the Wilson loop of 
an external qq-pair, at which the onset of string-breaking is normally expected. 
^^The correction ^ is very small: ii = -LZk£ ~ o.02. 

m ^ m ZTT m 



-15- 



string-breaking distance at temperatures close to the critical one and at low temperatures, 
smaller than ©(/jr). In both cases, the string-breaking distance increases with the increase 
of the temperature. We have also noticed that, in the case of spin-^ quarks, antiperiodic 
boundary conditions for fermions produce only corrections which fall off as the tail of the 
Gaussian distribution, as long as T < mT^/\/2. As a by-product, we have found in Ap- 
pendix B the temperature dependence of the string tension, which reproduces correctly 
the zero-temperature value. In Section 4, we have calculated (for scalar quarks) a cor- 
rection, generated by the second cumulant in the expansion (|4.1| ), i.e. by the dispersion 
of the chromoelectric field in the direction transverse to the string. This effect slightly 
diminishes the string-breaking distance. Using the Feynman variational method, we have 
then derived in Section 4 the leading quantum correction to this effect, produced by the 
deviation of the trajectory of a pair from the classical one. This effect further diminishes 
the string-breaking distance. 

The effects of dispersion of the chromoelectric field are small as long as the cumulant 
expansion is convergent, that is the case when the logarithm of the Landau-Ginzburg 
parameter is bounded from above according to ( [4.61) . This constraint is, however, always 
obeyed as long as at least the first pole in the Schwinger formula gives its contribution, 
that leads to an even more severe constraint ( |2.9| ). Although both constraints have been 
shown to relax at T ^ Tc, they do exist at T = 0. This fact necessitates to perform 
the average ([4.1| ) without the use of the cumulant expansion at all, that has been done 
in Section 6. As a result, the upper boundary on the logarithm of the Landau-Ginzburg 
parameter increases by a factor 4 [cf. eq. ( |6.1[ )]. Furthermore, a novel formula ( |6.3| ) for 
the rate of the pair production has been derived. At T = 0, some range of the Landau- 
Ginzburg parameter has been found, in which the string-breaking distance is larger than 
the one we had with the use of the bifocal approximation to the cumulant expansion in a 
factor varying between 2.84 and 9.10 Instead, at T — > Tc, this factor is smaller, namely 



it equals \/2 . Notice that such an analytic average of the exponent ( iA ) without the use 
of the cumulant expansion was only possible due to the explicit form of the Abrikosov- 
Nielsen-Olesen solution in the London limit. Apparently, studies away from this limit will 
require a numerical analysis. 

Within the approach, which treats pairs as holes in the confining pellicle, the new 
quantity on which the string-breaking distance is dependent is the constant entering the 
perimeter law of the Wilson loop of a produced pair. For typical values of the hadronic 
mass, string tension, and the vacuum correlation length, this dependence is, however, very 
weak. As for the dependence of the new expression for the string-breaking distance on the 
mass of a produced pair and on the string tension, it is the same as in the above-discussed 
case of the approach based on the Schwinger formula, where only the second cumulant is 
taken into account. Finally, the results for the string-breaking distance, obtained within 
all the above-mentioned approaches, are 0{N^) in the large-A'^c limit. 



^''This result indicates that, at least at these values of the Landau-Ginzburg parameter, the effect pro- 
duced by cumulants higher than the quadratic one is opposite and stronger than the result produced by 
the quadratic cumulant alone. 
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A. Some technical details 

Let us first present evaluation of the dispersion parameter ([4.5|). Integrating both sides 
of the equality ^Ko{mr) = f ^^J^^ ^1;+^^ o'^^r cPr we trivially have / d?rKQ{mr) = 

f ^ -a Hp) = In the same way 



frPK^( )=frfi MV'g e-(P+q)xx _ W £p _ 1 

jar o^mrj (2^)2 (p2 + ^2)(g2 + ^2) (^2 + ^2)2 2m^' ^ 

The parameter ([4.5|) then reads 



Next, we will discuss some details of derivation of Cvar- Taylor expanding Ko{z') in 
eq. ( ^.31) to the second order we have 

1 f j2^ ,2 / 



ml 



d'zd'z'e ^^"^v Ko{z)Koiz') 



27rT 



ml 



d'zK^iz)-^ I d^z^Kiiz) I d'z'{z' -z)^e 



-72 J I I 



-{z-z')2 



+ 



+ 



1m\ 



Sz 



'j^^^^^^K,{z) + '-^(Koiz) + K,{z)) 
z-^ Zz'^ 



Ko{z)x 



X d z'{z' - z)f^{z' - z)i,e 



-{z-z')2 



The second term on the r.h.s. of this equation apparently vanishes, while the third ones 
reads ^ / d^zKo [i(Ko + K2) - ^Ki] . Using the definition Cvar = ^ (^{£^) - {Ef) , we 
then arrive at the following intermediate result: 



c + 



%2l3 



00 

J dzKo [z{Ko + K2) - 2iCi] . 



Here, the addendum c, eq. (4^), is apparently produced by the term with no derivatives, 
while the other addendum is produced by the second-derivative term of the Taylor expan- 
sion. Finally, the integral over z can be calculated exactly. It reads 



-17- 



oo 

^ J (fzKl + J dzKo 



z[Ko + -^Ki]- 2Ki 



- f <fzKl = — 



where eq. (|A.lD on the last step has been used. This yields eq. (| 
We will finally present a proof of the formula 



oo 



jdx^,jdx'^G{{x-x'f) ^L- j diG{f), (A.2) 

C C 

1 

where L = J ds\x^{s)\ is the length of the contour C. Since the function G(x^) is rapidly 



decreasing [e.g. G{x^) = 2Di{0)Tg{\x\ + rg)e-l^'l/^f for the adopted Ansatz Dii^x"^) = 

L'i(0)e~l^l/"^s], we can Taylor expand x'^ = + t) as x^{s + t) ~ Xf^{s) + tXf^{s). In the 

1 1 

expression under study, f ds J dtx^{s)x^{s + t)G{t^xj^{s)), we therefore have x^{s)Xfj_{s + 



t) ~ i;^(s)+ti^(s)x^(s) = x^(s), where, at the last step, the proper-length parametrization 
xf^ = const has been fixed. Introducing instead of t the new integration variable ^ — t|x^|, 

1 

we have for the expression under study / ds\xi^{s)\ f d£^G{^'^). Finally, due to the rapid 



decrease of G(^^), the upper integration limit in the last integral can be replaced by infinity, 
that completes the proof of eq. ( |A.2|) . 



B. Temperature dependence of the string tension 

To find the dependence o"(T), notice that the string tension can be derived from the string 
representation of the partition function ||2^. An essential result of this representation is 
the following string effective action: 



2(7rt!)^y da^v{x) j da^,y{x')Dm^{x - x') 



where Dm{x) = mKi{m\x\)/{4Tr'^\x\) is the Yukawa propagator. The zero-temperature 
string tension then stems from this action according to the general (for this type of non- 
local string actions) formula [^: (Tq = (vrw)^ • J . This indeed equals 

^ mv/rriH 

the above-used value 2ttv'^L (following from the Landau-Ginzburg equations). At finite 
temperatures (smaller than the temperature of dimensional reduction), we rather have the 
following expression in terms of Matsubara frequencies: 



d^z y: 



my I mil 

To perform the integration, it is useful to transform the sum as follows: 
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+00 „ 

/ dtexp 



- - t (z2 + (my/3n)2) 



TJ^ ^ f dt 



my 



^2^ _ ]_+ (t^n/mv)" 



4t 



ttT 

mv\z\ 



+00 

E 



-|z|-y/l+(a;„/m-i/)^ 



Integration of this expression over d?z is now trivial and yields 

+00 



Performing now with the sum a transformation of the form ( |3.3| ) in the opposite direction, 
we obtain 

In particular, at T — > 0, one may approximate the sum by the zeroth term, that recovers 
the value ctq = 2ttv'^L. 
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